Highly anisotropic magnetic states of Co dimers bound to graphene-vacancies 
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The adsorption behavior and the magnetic states of cobalt atoms and dimers on single vacancies 
in a graphene sheet are investigated by means of relativist ic density functional calculations. It is 
found that local magnetic moments are formed in both cases, despite strong chemical binding. Of 
particular interest are kinetically stable isomers with two cobalt atoms attached to the same side of 
the graphene sheet. Magnetic bi-stability with an anisotropy barrier of about 50 meV is possible in 
this geometry. The feasibility of its preparation is discussed. 
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Nanoscopic magnets with large magnetic anisotropy 
(MA) are of interest both for fundamental research [T]-[3] 
and for technological applications [H [5]. Experimental 
activities to produce such systems were hitherto mainly 
focused on one of two known routes, the chemical syn- 
thesis of single-molecule magnets (SMM) or the physical 
deposition of magnetic atoms or small clusters on metallic 
substrates. Both routes resulted in systems with maxi- 
mum MA energies (MAE) somewhat below 10 meV per 
magnetic atom or per SMM: Single Co atoms deposited 
on a Pt(lll) surface show an MAE of 9 meV per Co 
atom pQ; an anisotropy barrier of 7 meV was reported 
for a complex of the Mn 6 family of SMM 0. 

Recently, the MA of four 3d transition metal dimers 
was predicted to be up to one order of magnitude larger 
than these figures [6j[7]. In the case of Co dimers, adsorp- 
tion on pristine graphene does obviously not harm their 
magnetic properties. Their MAE has been estimated by 
density-functional (DF) calculations [8] to be of the or- 
der of 100 meV. This finding adds a new flavor to the 
remarkable features of graphene |9j [10] (G) which has 
been praised [11] for its electronic, mechanical, optical, 
and thermal properties. 

Its structural robustness [11 , impermeability [12], 
chemical tolerance [13] . and reassemble properties [T4] 
distinguish G as a salient substrate material. As such, 
it has been applied to adsorb metal atoms [15] [16] and 
clusters [17]. Several theoretical studies confirmed that 
magnetic moments are to be expected on transition metal 
atoms and dimers including cobalt adsorbed on G [TSl - 
[2T] , Thus, C02-G seems to be suited for ultrahigh-density 
magnetic recording [5]. 

Unfortunately, the binding of transition metal atoms 
and dimers to pristine G is relatively weak. In the case of 
Co dimers, the DF binding energy amounts to about 0.7 
eV. [8] However, a yet smaller activation energy of 0.24 
eV (upper limit) for surface migration was measured for 
gold atoms on graphite [22] and migration barriers in the 
range of 0.01 to 0.8 eV were evaluated [T51 [23] bv DF 
calculations for transition metal adatoms on G. At room 
temperature, Co atoms are hence expected to float on the 



carbon layer [23 . Indeed, STM images of Co atoms on 
pristine G were only obtained at 4.2 Kelvin [24], hitherto. 

Defect structures of carbon are highly reactive. Thus, 
one possibility to bind transition metal atoms or dimers 
more strongly consists in a deliberate creation of struc- 
tural defects, e.g. single vacancies with dangling bonds, 
in the G sheet before the adsorbate is deposited. An 
established experimental route to produce isolated va- 
cancies in graphene-like systems is irradiation with low- 
energy Ar + ions [25]. After annealing, most defects pro- 
duced in this way were found to be single vacancies [26] . 

For the case of Co atoms, a binding energy to a sin- 
gle vacancy in graphene (SVG) of about 8 eV was pre- 
dicted [23], by far large enough to ensure stability at 
room temperature. Strong chemical binding, on the other 
hand, often spoils magnetism. In particular, orbital mag- 
netism which is chiefly responsible for the MA is known 
to be quenched by hybridization and crystal fields. While 
the existence of spin magnetism was predicted for atomic 
impurities of several transition metals in G including Co- 
SVG [23j [27] , we are not aware of related investigations 
concerning orbital magnetic properties. It is thus both 
demanding and interesting to find magnetic graphene- 
based nanostructures that are structurally stable on the 
one hand and show a large orbital moment or even a large 
MA on the other hand. 

Here, we present C02-SVG as an example of a nanos- 
tructure that possibly fulfills both contradicting require- 
ments. The related system Co-SVG is demonstrated to 
be a counterexample where the strong binding suppresses 
the MA. As a specific intricacy, the desired large MA 
of C02-SVG is not present in the ground state but in 
kinetically stable isomers which can be prepared in the 
most simple way by atom- wise deposition. The suggested 
preparation can be considered as a combination of the 
two known routes to produce nanomagnets with large 
MA, the chemical (SMM) and the physical (adsorbate 
deposition) one. 

The DF calculations were performed with the all- 
electron full-potential local-orbital (FPLO) code [28] . 
version 9.00-34. The exchange-correlation energy func- 
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TABLE II: Orbital-resolved occupation numbers of 3d spin- 
down states of r-Co for the three electronic isomers of the 
geometry (A) and for geometry (C). The quantum numbers 
m refer to real spherical harmonics. 



FIG. 1: (Color online) Four isomers (A)-(D) of Co 2 -SVG ob- 
tained by relaxation of five initial geometries [30] . Large blue 
spheres: cobalt atoms; small gray spheres: carbon atoms. 
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TABLE I: Relative total energies, .Etot, and total spin mo- 
ments, /i s , of the structural isomers of C02-SVG, Fig.[l] E sp i n 
denotes the stabilization energy of spin magnetism, i.e., the 
energy difference between non-magnetic and magnetic state. 
Structure (A) can be stabilized in three different electronic 
states: (AO), (AI), and (All). 



tional was evaluated using the standard parameteri- 
zation [29] of the generalized gradient approximation 
(GGA). Orbital magnetic properties were calculated in 
a four-component fully relativistic mode. To obtain an 
upper estimate for the MAE [8] , orbital polarization cor- 
rections [31] (OPC) were applied in some of the calcula- 
tions. The technical parameters of the calculations are 
described in the Online Supplement [30] . 

Recent DF results indicate that single Co atoms on 
pristine G are magnetically bistable with an MAE-barrier 
in the 10-meV range [8]. Further, using GGA it has al- 
ready been found that Co atoms bound to an SVG keep a 
spin moment of 1 /iB , while Fe and Ni atoms become non- 
magnetic [23] [27] . These predictions tempted us to check 
the magnetic anisotropy of Co- SVG, before considering 
the more complex case of Co dimers. The values of spin 
moment, /i s (Co-SVG) = 0.92 /iB, energy needed for Co 
detachment, £?co-det(Co-SVG) = 8.0 eV, and Co-C bond 
length, d c °- c (Co-SVG) = 1.76 A, found in our calcula- 
tions [30] confirm the available literature data [23 , 27, 32 . 
By checking the orbital magnetic properties, we how- 
ever find negligible magnitudes of the orbital moment, 
/iz(Co-SVG) < 0.01 /i B , and of the MAE (< 0.1 meV) in 
both GGA and GGA+OPC approaches. Thus, we con- 
clude that the orbital magnetism of single Co atoms is 
quenched by binding to an SVG. 

There are several structures for the interaction of two 
Co atoms with SVG imaginable. We considered five 



initial geometries [30] . Relaxation of these geometries 
resulted in four (meta)stable structures, (A) to (D) in 
Fig. [l] In the case of (A), we found three different 
electronic isomers (AO), (AI), and (All). They are dis- 
tinguished by specific 3d-orbital occupations at the "re- 
mote" Co atom (r-Co) which is not directly bound to the 
vacancy. The relative total energies and the spin mo- 
ments of all isomers are given in Tab. [TJ their structural 
details are compiled in the Online Supplement [30] . 

All identified C02-SVG isomers are magnetic with spin 
S = 1. This is a notable difference to the S = 2 states of 
free Co dimers [6] and of C02-G [8]. The lower spin state 
of C02-SVG mirrors the stronger chemical bonding of Co 
at the vacancy compared to the pristine carbon hexagon. 

Table [il] compiles the 3d spin-down orbital occupation 
numbers of r-Co for (A0-II). In all three isomers, the 
spin moment resides on r-Co and the Co atom at the 
vacancy (v-Co) is nearly spin-compensated. Thus, the 
spin-up channel of r-Co is almost completely occupied. 
State (AO) is lower in energy than (AI) and (All), but 
its charge density lacks axial symmetry. Using a GGA+U 
mode with U, J — >• we were able to stabilize all three 
states in a scalar relativistic GGA calculation for a free 
Co atom. The energy sequence is now (AO, free) - (All, 
free) at 0.1 eV - (AI, free) at 0.2 eV. The correct atomic 
ground state has an axially symmetric charge density like 
(AI, free) or (All, free), but unlike (AO, free). If spin- 
orbit interaction (s-o) is switched on, we find L = for 
(AO, free), L = 3 for (AI, free), and L = 2 for (All, 
free). Considering Hund's second rule, the best GGA 
approximation to the atomic ground state is (AI, free), 
the second best is (All, free). Both facts, its broken 
axial symmetry and its incompatibility with Hund's sec- 
ond rule, lead us to disregared the state (AO, free) as 
unphysical. Its low energy is due to self-interaction of 
the non-spherical part of the charge density, as discussed 
earlier by Brooks et al. [34 . Bonding of the free Co 
atom to Co-SVG slightly shifts the energies of the three 
electronic isomers, but hardly changes their orbital oc- 
cupations. We therefore exclude (AO) from the further 
discussion. 

Structure (B), where the two Co atoms are located on 
opposite sides of the G plane, is the ground-state iso- 
mer. However, we will now provide arguments that it is 
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FIG. 2: (Color online) Diffusion barrier between C02- 
SVG(AI) and Co 2 -SVG(B). The inset shows the considered 
structure. The height of the Co atom next to the vacancy, 
h Go , was fixed (abscissa of the figure). Black (open) symbols: 
space group 156 (P3ml), red (full) symbols: space group 1 
(PI). Relaxation was allowed for the height of the three 
neighboring carbon atoms, h c of Panel (a), their lateral posi- 
tion, and the distance between the two Co atoms (not shown). 
Panel (b) displays the total energy relative to the energy of 
Co2-SVG(AI). The two minima in the energy curves are re- 
lated to the two isomers (AI) and (B), the diffusion barrier 
amounts to 2 ... 2.5 eV, estimated by extrapolation of the left 
branch. 



technically easier to produce a structure with both Co 
atoms on the same side of the carbon layer and that this 
situation is kinetically stable. 

Deposition of small Co amounts on a graphene-like sur- 
face with aforehand-prepared single vacancies will yield 
Co-SVG under ambient conditions, since Co atoms are 
mobile on pristine surface areas. A fraction of the Co 
atoms may agglomerate [Co-G —> G + hep-Co], but the 
related energy gain per Co atom [30 j, Ec - ag g = —4.4 
eV, has a smaller absolute value than the energy gained 
by binding to the SVG [Co-G + SVG -> G + Co-SVG], 
£7 b (Co-SVG) = -6.8 eV. We further note, that Co- 
SVG is expected to be stable against migration at room 
temperature. Related activation energies for Pt-SVG 
and Au-SVG were found to be close to 2.5 eV by high- 
temperature transmission electron microscopy [33] . 

If Co-deposition is continued after saturation of the 
vacancies, a large fraction of the mobile Co atoms will 
attach to Co-SVG and form Co2-SVG(AI) or one of 
its quasi-degenerate isomers (All) or (C). The reaction 
[Co-G + Co-SVG -> G + Co 2 -SVG(AI)] gains an en- 
ergy [30^ of E h (Co 2 - SVG) = -1.2 eV. Agglomeration 
of the Co atoms with a gain of -4.4 eV would be ener- 
getically yet more favorable, but obviously the formation 
of C02-SVG is kinetically more likely at low deposition 
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TABLE III: Results from fully relativistic GGA calculations 
for Co 2 -SVG(A) and Co 2 -SVG(C). The notations 001 and 
100 refer to magnetic moment orientations perpendicular and 
parallel to the C plane, respectively. Magnetic spin moments 
/i s and orbital moments fii are projected on the cobalt atoms 
next to the vacancy (v-Co) and on the other Co atom (r-Co) 
and given in jib- The MAE is the difference between total 
energies obtained for the two chosen moment orientations, 
MAE = £7(100) - £7(001), and refers to the whole structural 
unit with two Co atoms. The in-plane anisotropy is small and 
can be neglected. Results in parentheses refer to GGA+OPC 
calculations. 



rates and sufficient supply of Co-SVG sites. The sub- 
strate temperature should be low enough to avoid de- 
tachment of Co from C02-SVG (below 300 K, estimated 
from |£ b (Co 2 - SVG) I = 1.2 eV), but high enough to 
ensure mobility of the deposited Co atoms on G. 

Finally, we check if Co2-SVG(AI) is kinetically sta- 
ble against diffusion of the Co atom next to the vacancy 
to the other side of the G sheet [Co 2 -SVG(AI) -> Co 2 - 
SVG(B)]. Fig. [2] presents the geometry and the energy of 
structures along the diffusion path. We find an energy 
barrier of more than 2 eV, large enough to ensure stabil- 
ity beyond detachment of the upper Co atom. This large 
barrier arises from the fact that the Co atomic radius is 
much larger than the radius of the missing carbon atom. 
Thus, Pauli repulsion on the one hand and stiffness of 
the carbon plane on the other hand hinder the diffusion. 
We conclude that preparation of C02-SVG with both Co 
atoms located at the same side of the C plane is possible 
and that this structure is stable up to about 300 K. 

Considering the size of the discussed self-energy error 
for systems containing almost free atoms, we are not able 
to judge if (AI), (All), or (C) is the ground state isomer. 
We however tend to exclude geometry (D), which lies 
0.4 eV above (AI), from the further discussion. We also 
note that the magnetic state of (AI), (All), and (C) is 
more stable than that of the other isomers, cf. the spin 
stabilization energies in Tab. [I| 

We now turn the attention to orbital magnetic prop- 



erties. Tab. Ill compiles results for the kinetically stable 
low- lying isomers. In the case of (AI), any self-consistent 
GGA+OPC calculation converged to the state (All). 
This might be an indication that (AI) is not realized in 
nature. Most remarkably, the MAE of (All) and of (C) 
in the GGA+OPC approach is predicted to be similarly 
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high as that of C02-G [8], though in the present case the 
Co dimer is much more tightly bound to the substrate 
and thus has a much better chance to be accessible to 
experimental realization. 

Strong chemical bonding claims its price, though. The 
atomic-site projected magnetic moments of the Co atom 
next to the vacancy are almost completely quenched, 
Tab. Ill However, large spin moments (S = 1) persist 
on r-Co due to the low coordination of this atom. For 
(All) and for (C) in OPC, also the orbital moment is 
large, L = 2. Magnetic moment orientation perpendicu- 
lar to the easy (001) axis quenches the orbital moment 
but requires a large energy, the MAE, which is related to 
the orbital moment reduction [35j [36] . 

Among the two Co atoms, a distribution of tasks takes 
place: v-Co atom is strongly bound to SVG and serves 
as an anchor for r-Co; the remote r-Co, on the other 
hand, preserves atomic-like properties and carries those 
electronic states which produce an ultimately high MAE. 
From this distribution of tasks it is clear that atoms of el- 
ements other than cobalt could be used as anchor atoms 
and/or remote atoms. Thus, C02-SVG is most likely 
only one example of a whole class of unknown highly 
anisotropic nanomagnets. 

To summarize, we suggest that deposition of Co atoms 
on SVG can be used to prepare a structure with two Co 
atoms attached to each vacancy. We demonstrated its 
kinetic stability and found strong indications that it will 
show an exceptionally high MA. Single Co atoms on Pt 
also show large MA pQ , but their technological use might 
be restricted by a high spin tunneling probability [3 . For 
the present system, the hybridization of the magnetically 
active atom with its surrounding is much smaller than 
in the case of a metallic substrate, since the chemical 
bonding is mediated by a single anchor atom. Thus, C02- 
SVG might even be of technological interest for magnetic 
data storage. 

We finally note that C02-SVG should be considered as 
one didactic example of a whole class of possible struc- 
tures, consisting of an organic substrate with dangling 
bonds and a few metal atoms. The "remote" atoms show 
strong orbital magnetism, while the "anchor" atoms me- 
diate the bonding to the substrate. The production of 
such systems could make use of a self-assembling step: 
the saturation of the dangling bonds by anchor atoms. 
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Supplementary material 

This Supplement reports technical details of the cal- 
culations, results of the optimization of geometry and/or 
spin state, definitions of reaction energies and conver- 
gence checks for the simulation cell size. 



Technical details of the calculations 

The convergence with respect to the /c-space grid den- 
sity was carefully checked. Using a linear tetrahedron 
method with Blochl corrections, the final calculations 
were performed with 9x9x1 /c-points in the full Bril- 
louin zone (geometry optimization) or with 30 x 30 x 1 
/c-points (orbital magnetism). The basis set comprised 
Co (3s, 3p, 3d, 4s, 4p, 4d, 5s) and C (Is, 2s, 2p, 3s, 3p, 
3d) states. The modified G sheet was modeled with a 
3x3 supercell of fixed edge length 7.392 A in the x-y 
plane and with 12 A spacing along the z direction. Cal- 
culations with a 4 x 4 cell and with 16 A spacing did not 
show any significant difference, see below. 

Geometry optimization was carried out using forces 
evaluated in a scalar relativistic mode with a convergence 
criterion of 1 me V/A. Fixed spin-moment calculations 
were employed to exclude the existence of spin states not 
found by the unrestricted self-consistent calculations. 

Optimization of geometry and/or spin state 

The geometry and/or the spin magnetic state of eight 
different systems were optimized with respect to the total 
GGA energy: a free Co atom, a free Co dimer (C02), bulk 
Co metal in the hexagonal close-packed structure (Co- 
hcp), a Co atom on pristine graphene (Co-G), a Co dimer 
on pristine graphene (C02-G), graphene with a single va- 
cancy (SVG), Co-SVG, and Co 2 -SVG with four identified 
structural isomers, called (A, B, C, D), see Fig. [ST] and 
Fig. [l] of the main text. In the isomer C02-SVG (B), 
the two Co atoms are located on opposite sides of the 
graphene sheet. Thus, this isomer is also referred to as 
Co-SVG-Co. 

All graphene-derived systems were treated in a 3 x 3 
structural unit cell. For comparison, the total energy of 
pristine graphene (G) is computed in a fixed geometry, 
according to the fixed unit cell size used in the other 
calculations. All energies and spin moments in Tabs. [Si] 



and SII refer to the considered structural unit, i.e., to 18 



carbon atoms for G, to 18 carbon atoms and 1 (2) cobalt 
atom(s) for Co-G (C02-G), to 17 carbon atoms for SVG, 
and to 17 carbon atoms and 1 (2) cobalt atom(s) for 
Co-SVG (Co 2 -SVG). 

All calculations were performed with the code and set- 
tings defined in the main text, exceptions are mentioned 
in the column 'remarks'. For Co atoms and dimers, the 




FIG. SI: Four isomers (A)-(D) of Co 2 -SVG obtained by 
relaxation of five initial geometries, see Fig. |S2] Blue spheres 
denote cobalt atoms and gray spheres denote carbon atoms. 
Fig.[I]of the main text shows the same structures in a different 
view. 



molecular mode with free boundary conditions was ap- 
plied. All other calculations were performed with three- 
dimensional periodic boundary conditions. 

In the atomic calculations, the 3d and 4s states turn 
out to be degenerate. The lowest GGA energy corre- 
sponds to a configuration 3d 7,6 4s 1,4 . Spin-orbit coupling 
has a stronger influence on the atomic total energies than 
on the total energies of atoms in a cluster, molecule, 
or bulk compound. We included related data obtained 
in fully relativistic mode for the Co atom and the Co 
dimer in Tab. [Sjbut used the values obtained in scalar 
relativistic mode to calculate binding energies for the 
sake of comparability with other published (scalar rel- 
ativistic) DF energies [23 . Binding energies related to 
atomic Co would have to be corrected by the amount 
E{Co - atom) - E r {Co - atom) = 0.08 eV. The correc- 
tion related to a Co dimer amounts to less than 0.02 eV, 
see Tab. ED 

Fig. [S2] shows the five chosen initial structures used 
to identify isomers of C02-SVG by structural relaxation. 
The initial SVG structure was always prepared by remov- 
ing one carbon atom from the pristine graphene sheet 
without relaxation. The two Co atoms were initially 
placed in the following positions to form the initial struc- 
tures (I) to (V): 

(I) both above the vacancy site at heights of 1.5 and 3.5 
A; 

(II) above and below the vacancy site at heights of 1.5 
and -1.5 A; 
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performed in space group 1, the other two, (AI) and 
(All), in space group 156, using different starting con- 
ditions. 





III 



IV 



FIG. S2: Five chosen initial geometries (I)-(V) of Co 2 -SVG. 



(III) above the centers of two pentagons next to the va- 
cancy at the same height of 1.5 A; 

(IV) above the vacancy site and above the center of a 
neighboring hexagon at the same height of 1.5 A; 

(V) above the vacancy site and above the center of a 
pentagon next to the vacancy at the same height of 1.5 

A. 

The initial geometries for Co-G and Co- SVG were pre- 
pared by positioning the cobalt atom above the center of 
a hexagon at a height of 1.5 A and above the vacancy site 
at a height of 1.5 A, respectively. In the case of C02-G, 
the Co atoms were placed above the center of a hexagon 
at heights of 1.5 and 3.5 A. 

A summary of results for all optimized geometries is 



given in Tabs. [Si] and |SII| The calculated atomic dis- 



tances d and heights h in Tabs. SI and SII refer to nearest 
neighbor distances and to the distances to the x-y plane, 
respectively. In particular, d c ~ c denote the distances 
between pairs of the three carbon atoms next to the va- 
cancy, h c is the displacement of these atoms from the 
x-y plane, d Co ~ c are the distances between the cobalt 
atom next to the vacancy and the related carbon atoms, 
and h Co is the distance of the cobalt atom next to the 
vacancy from the x-y plane. 

A notable structural deformation was found at the va- 
cancy site for the Co2-SVG-isomers (A), (C), and (D), 
where the three carbon atoms next to the vacancy are 
pulled out up to 0.6 A from the graphene plane. 

In the case of structure (A) we found three electronic 
isomers. Isomer (AO) was obtained if the calculation was 



Reaction energies 

We define reaction energies of the following processes, 
referred to in the main text: 

Energy needed for the detachment of Co atoms from SVG 
(reaction Co-SVG Co-atom + SVG), 

^Co-det(Co-SVG) = 

£(Co-atom) + £(SVG)- 

-£(Co-SVG) = 7.96 eV. (1) 

Energy gained by agglomeration of Co atoms on graphene 
(reaction Co-G — » G + hep-Co), 

^Co-agg = 

£(G) + £(Co-hcp)- 

-£(Co-G) = -4.43 eV. (2) 

Energy gained by binding of Co atoms on graphene to 
SVG (reaction Co-G + SVG -> G + Co-SVG), 



£ b (Co-SVG) 
£(G) + £(Co-SVG)- 
-E(Co - G) - £(SVG) 



-6.76 eV . 



(3) 



Energy gained by binding of Co atoms on graphene to 
Co-SVG. For reasons explained in the main text, we take 
isomer (AI) as reference (reaction Co-G + Co-SVG — » G 
+ Co 2 -SVG(AI)), 

£ b (Co 2 -SVG) = 
E(G) + E(Co 2 - SVG(AI))- 
-E(Co - G) - E(Co - SVG) = -1.15 eV. (4) 

Convergence of the simulation cell 

We checked the magnetic moment and binding energy 
of Co-SVG with a 4 x 4 supercell and found that the 
ground state has a magnetic moment of 0.99 //b- The Co- 
detachment energy of this structure amounts to 7.99 eV. 
These values are almost the same as those obtained with 
a 3 x 3 supercell, 0.92 /iB and 7.96 eV. 

In the case of SVG, the 3x3 supercell yields equal 
distances between the carbon atoms next to the vacancy 
and a vacancy formation energy, 

£vac(3x3) = 

£(SVG,3 x 3) - £(G,3 x 3) • 17/18 = 7.95 eV . (5) 



If a 4 x 4 supercell is used, the so-called 5-9 structure 
is obtained where two of the carbon atoms next to the 
vacancy form a loose pair. Accordingly, the vacancy for- 
mation energy is about 0.1 eV smaller than in the former 

case, 

£vac(4x4) = 

£(SVG,4 x 4) - £(G,4 x 4) • 31/32 = 7.84 eV . (6) 

The dependence of £^(SVG) on the size of the simu- 
lation cell is compensated by an opposite dependence of 
E(Co - SVG). Thus, £ Co -det(Co - SVG) is almost in- 
dependent of the cell size. We conclude that the reaction 
energies are reliable at a level of 0.1 eV, if 3 x 3 supercells 
are used. 
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system 


space 
group 


geometry 
parameters 

[A] 


total 

energy 

[eV] 


spin 

moment 

[Mb] 


remarks 


Co atom 
Co atom 






E = -37916.34 
E T = -37916.42 


3.00 
3.00 


scalar relativistic 
fully relativistic 


Co 2 
Co 2 




d c °- Co = 1.989 
^co-co = L9g9 


E = -75836.44 
E r = -75836.45 


4.03 
4.03 


scalar relativistic 
fully relativistic, 
distance fixed 


Co-hcp 


194 

(P63/mmc) 


a = 2.5071 
c = 4.0695 


E = -37921.97 


1.62 


experimental geometry 
/c-mesh 96 x 96 x 60 
scalar relativistic 


G 


156 (P3ml) 


dg~ c = 1.4226 


E = -18663.46 





experimental geometry 
/c-mesh 60 x 60 x 1 
scalar relativistic 



TABLE SI: Compilation of reference energies and other prop- 
erties of Co atoms, Co dimers, bulk Co, and pristine graphene. 



system 


space group 


geometry parameters 

[A] 


total energy 
[eV] 


spin 

moment 

[Mb] 


Co-G 


183 (P6mm) 


d Co-c = 2 in 
h Co = 1.560 


E = -56581.00 


1.06 


C02-G 


183 (P6mm) 


d c °- Co = 2.054 
d Co-c = 2250 

h Co = 1.743 


E = -94500.57 


4.06 


SVG 


1 (PI) 


d c ~ c = 2.545, 2.546, 2.545 
h c = 


E = -17618.65 


1.04 


Co-SVG 


1 (PI) 


d Co-c = x 76 2 5 x 76 2 5 L762 
h Co = 1.144 


E = -55542.96 


0.92 


C02-SVG 
(AO) 


1 (PI) 


d c °- Co = 2.250 

d Co-c = x x x 763 

/i Co = 1.195 

/i c = 0.270, 0.275, 0.275 


E = -93461.76 


2.00 


C02-SVG 
(AI) 


156(P3ml) 


d c °- Co = 2.280 
d Co-c = 1764 

h Co = 1.200 
h c = 0.275 


E = -93461.65 


2.00 


C02-SVG 
(All) 


156(P3ml) 


d c °- Co = 2.246 
d Co-c = 1764 

h Co = 1.200 
/i c = 0.276 


E = -93461.48 


2.00 


C02-SVG 
(B) 


1 (PI) 


d Co-c = l g76 ^ l g76 ^ 1S76 

h Co = 1.171 

h c = 0.0, 0.0, 0.0 


E = -93462.68 


1.71 


C02-SVG 
(C) 


1 (PI) 


d c °- Co = 2.303 

d Co-c = 176Q ^ l 7m ^ i m 

h Co = 1.238, 2.566 

h c = 0.232, 0.231, 0.586 


E = -93461.54 


2.00 


C02-SVG 
(D) 


1 (PI) 


d c °- Co =2.338 

d Co-c = x 76 q 5 1>784> x 763 

h Co = 1.240, 1.986 

/i c = 0.154, 0.312, 0.445 


E = -93461.27 


1.75 



TABLE SII: Compilation of geometries, total energies, and 
spin moments of Co-G, Co 2 -G, SVG, Co-SVG, and Co 2 -SVG 
structures, obtained in scalar relativistic mode. 



